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The proximity effect and the Josephson current in a superconductor-insulator-ferromagnet- 
superconductor (SIFS) junction are investigated within the framework of the quasiclassical Eilen- 
berger equations. This investigation allows us to compare the dirty and the clean limits, to in- 
vestigate an arbitrary impurity scattering, and to determine the applicability limits of the Usadel 
equations for such structures. The role of different types of the FS interface is analyzed. It is shown 
that the decay length and the spatial oscillation period 2-k£,2 of the Eilenberger function may 
exhibit a nonmonotonic dependence on the properties of the ferromagnetic layer such as exchange 
field or electron mean free path. The results of our calculations are applied to the interpretation of 
experimentally observed dependencies of the critical current density on the ferromagnet thickness 
in Josephson junctions containing a Ni layer with an arbitrary scattering. 



PACS numbers: 74.45.+C, 74.50.+r, 74.78.Fk 

I. INTRODUCTION 

Superconductor- ferromagnet-superconductor (SFS) 
Josephson junctions are a subject of intensive theoretical 
and experimental studies^"— . In particular, the question 
of the applicability range of predictions from dirty and 
clean limit theories and the treatment of the crossover 
between these two limits has been recognized as an 
important problem for the theoretical description of SFS 
structures^. 

For the majority of experimental realizations of SFS 
structures'^— the exchange energy H of the ferromag- 
netic materials is rather largeiSiii. As a consequence, 
the characteristic magnetic length = hvf/2H < 
where ^/ is the electron mean free path and w/ is the 
Fermi velocity in the F layer—"—. Under this condition 
the numerous theoretical predictionsi"— based on the Us- 
adel equationsi' have a rather restricted range of valid- 
ity; i.e. a more general approach based on the Eilenberger 
equationsii has to be developed. 

A simple expression for the critical current density Jq, 
which is valid in the clean limit, was derived in Ref. llSl . 
Still the analysis of the Eilenberger equations for more 
general cases, remained a difficult problem. A significant 
progress along this direction has been achieved a decade 
ago^, where the solution of the Eilenberger equations for 
arbitrary scattering has been expressed in an integral 
form^. It was supposed^ that the SF interface trans- 
parency D is small enough, providing the opportunity to 
use the linearized equations. Thus, a general expression 
for the Josephson junction supercurrent has been derived 
and used for numerical evaluations. The integral rep- 
resentation of the solution of the Eilenberger equations 
also permits to reproduce analytical expressions for Jc 
obtained earlier within both clean and dirty limits' d^i^^ . 
Recently, the same results have been achieved using the 



Ricatti parametrization of the one-dimensional Eilen- 
berger equations^. However, the obtained expressions^ 
are so complicated that they arc difficult analyze and 
use in practice. Although the use of a one-dimensional 
equation' significantly simplifies the problem, it contains 
a non-controllable assumption about the insignificance of 
the angular distribution of the Eilenberger function. 

Despite of this progress, it is still not clear within which 
range of parameters and with what accuracy it is possible 
to use simple expressions for the clean'' and dirtj*i' lim- 
its. The question about the influence of transport prop- 
erties of SF interfaces is also still open. Since the answer 
to these questions is rather important for experimental- 
ists, we formulate here a problem for a particular case 
of superconductor- insulator- ferromagnet-superconductor 
(SIFS) junctions. In comparison with SFS, in SIFS struc- 
tures superconductivity is induced in the F layer only 
from one S electrode. This essentially simplifies the anal- 
ysis compared to that for SFS junctions. 

On the one hand, SIFS junctions are interesting by 
themselves as potential elements for superconducting 
logic circuits^'"— because it is possible to vary their crit- 
ical current density in the tt state within a rather wide 
range, still keeping a high JcRn produc t^^'^'^ , where Rn 
is the junction normal resistance per square. Further- 
more, in comparison with SFS junctions, SIFS junctions 
have very small damping which decreases exponentially 
at T — > 0. This makes them useful for superconducting 
circuits where low damping is required. On the other 
hand, an SIFS junction represents a convenient model 
system for a comparative study of O-tt transitions for an 
arbitrary ratio of characteristic lengths in the F layer: the 
F-layer thickness df, the mean free path £f, the charac- 
teristic magnetic length , and the nonmagnetic coher- 
ence length ^0 ~ hv f /2'KkBT, where T is the tempera- 
ture. 
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To calculate Jc of an SIFS junction, it is sufficient to 
study the proximity effect in the FS bilayer and to cal- 
culate the magnitude of the Eilenberger functions at the 
IF interface. For simplicity we will restrict ourselves to 
the study of the situation when the anomalous Green's 
function induced in the F layer is small enough, permit- 
ting us to use the linearized Eilenberger equation for the 
description of the superconducting properties induced in 
the fcrromagnet. Such an approximation is valid if the 
FS interface has a small transparency or if T is close to 
the critical temperature of the S electrodes. 

This article is organized as follows. In sec. [Tl] 
we describe our model based on the linearized Eilen- 
berger equation supplemented with Zaitsev boundary 
conditions24 for an SIFS junction. Different types of FS 
boundaries are analyzed. Section [IIII presents dependen- 
cies of anomalous Green's functions and the critical cur- 
rent density on the F-layer parameters. The comparison 
with experimental results is presented in Section llVl Sec- 
tion |V] concludes this work. The calculation details can 
be found in the Appendix. 



sign(i:j) d ^ 



(7) 



Here oj = TiT{2n + 1) are Matsubara frequencies, 9s j 
are the angles between the FS interface normal and the 
direction of Fermi velocities Vg.f in the S and F layers, re- 
spectively, /i/.s = cos 9 f^s, A is the superconducting order 
parameter, which is assumed to be zero in the fcrromag- 
net, Tsj = Isj/vsj are electron scattering times, and 
corresponding wave vectors ks = 2|a;|Ts-|-l, kf = 2[|cj| + 
iiJsign(w)]r/ + 1 = 2\u\tf /vf +i2Hif /vj sign(a;) + 1 = 
1 + ^f/^u + is,igTL{uj)lf /^H, = Vf/2 \uj\ . Wc use the 
units where h = 1 and fc^ = 1- Here wc also take into 
account that the normal Eilenberger function in the fcrro- 
magnet g « sign(a;). We also neglect multiple reflections 
from the FS and the IF interfaces, which is reasonable if 
df > if, i.e. the ballistic regime is not considered. The 
fully ballistic case with a double barrier SIFIS junction 
was examined in Ref. (25| . 

Equations (I3|)-(l7l) must be supplemented by the 
boundary conditions. At a; — ^ oo the function should 
approach its bulk value 



II. MODEL 

A. Proximity effect in the FS electrode 

To analyze the proximity effect in an FS bilayer for 
arbitrary values of the electron mean free paths is and 
if in the S and F layer, respectively, it is convenient to 
introduce the following functions 

^+ = ^[f{x,es,u) + f+{x,9s,oj)], (1) 



-i>-^-[fix,9s,u;)-f+ix,9s,u;)], (2) 

where f(x,9s,uj) and {x,9s,uj) are the quasiclassical 
Eilenberger functionSiii. Then we rewrite the Eilenberger 
equations for the S layer, located at < x < oo, in the 
form: 



kl 92 2Ar, + ($+) 

— $j_ — <I>^ = ^ — — 

i^sf^^dx^ + + 2\uj\ts + 1 ' 



where 



($+) = f 
Ja 



is signjuj) 

ks dx ^ 



(3) 



(4) 



(5) 



The insulating layer of negligible thickness is located at 
X — ~df. In the F layer (located at ~df < x < 0) the 
equations have the form 



«2 ,,2 Q^2 ^+ 



i'j^'j dx^ 



(£±) 

kf ' 



(6) 



(8) 



where Ao is the magnitude of superconductor order pa- 
rameter far from the SF interface. At x - 
boundary condition 



df the 



dx 



'f+{-df,9f,u)^0 



(9) 



guarantees the absence of a current across the insulating 
layer. 

The boundary conditions at the FS interface strongly 
depend on its transport properties. Below we will assume 
for the FS interface transparency -D(/i/) ^ 1. Then to 
the first approximation we may neglect the suppression 
of superconductivity in the S electrode and rewrite the 
Zaitsev boundary conditions^ in the form 



Ao 



if d 

kf dx vAf+^ 



(10) 



From the structure of Eilenberger equations (|3])-(l7]) 
and the boundary conditions ([8|)-(fT0|) it follows that 



(11a) 
(lib) 



<^+{9s,L0) = '^+{9s,-Uj), 

<^+{9f,^) = <^X{9f,-u), 



^^{9s,io) = -^-{9s,-oj), 
$_(0/,c^) = -$_(0/,-c^). 



(12a) 
(12b) 



These symmetry relations permit us to consider the so- 
lution of the boundary problem formulated above only 
for Lu > Q. 
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B. Solution of the Eilenberger equations in the F 
layer. 

It is convenient to look for a solution of Eilenberger 
equations in the F layer ([7]) in the form 

(13) 

which automatically satisfies the boundary condition ^ . 
The relation between the coefficients Qm, and B can be 



found by substitution of the ansatz (|T3)) into the Eilen- 
berger equation ([5]). Multiplying the obtained equations 
by cos[7rfc(a; + df)/df], and integrating them over a;, one 
can easily find the relation between the coefficients Q™, 
and B, see the Appendix. B = B{iJLf) can be found 
from the boundary condition (1101) . Thus we arrive to the 
following expression for the Eilenberger function in the 
fcrromagnet 



Ao 



E 



Trm{x + df) 



sinh 



cosh q 



X + df 



(14) 



where 



dfkf 



1 1 1 



Then the Eilenberger function averaged over the angle 9f has the form 



Ao 



E 



•q A/(m) 



\kf - ^arctanf^ 

I / Trm \ q 



COS 



mn{x + df) 



(15) 



There are several interface models with a different D{iif ) dependence. 

First, the FS interface may be represented by a thick diffusive barrier. Then the incident electron scatters in any 
arbitrary direction with equal probability, independent on the incident angle, i.e.. 



D{iif)^D,. 



(16) 



The second model considers an FS interface with a potential barrier, which appears due to different Fermi velocities 
in the S and F layer. The transmission coefficient in this case has the form 



AVsHsVfUf 4:Vf 



fif = D2fif, if Vs > vf 



{Vsfis+Vf^if) 

It is necessary to note, that the incidence angle and the reflection angle are related by the expression 

Vs sin 9s = Vf sin 9 / . 



(17) 



(18) 



Therefore, only the electrons which are at almost normal incidence to the interface (fj,s ~ 1) may penetrate through 
the barrier. This model seems to be mostly reasonable for description of recent experiments^ with Nb electrodes 
{vs ^ 5 . . .6 ■ lO^cm/s) and 3d ferromagnets or its Cu alloys^l (for Ni Vf « 2.8 ■ lO^cm/s). 

Third, one can model the FS interface as a high and narrow (J-function like) potential barrier between two metals 
with close Fermi-velocities. In this case 



AVsfJ.sVffJ-f 



{Vsfls + Vfflf + W) 



4:VsVf 2 



D 



if 



Vf, 



(19) 



where W is the strength of ^-function like barrier. 

For the above mentioned three models of the FS interface we have three possible expressions for with different 
average values ($+), sec (|A7|) - (|A9|) for details. 
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First, for D{nf) = Di 

AoDi 



E 



COS 



fcf - -2- arctan f ^l™") 



Second, for D{p,f) = /i/-D2 



Third, for D{fif)^ i^ijDs 



^ (^1)" fl - ^ arctan^) cos f 



■ arctan 



E 



(-ir 1- 



■? :j In 



( ■Km \ 
1 J 

(^ + l))cos(^ 



-2- arctan f 

7rm \ ^ J 



sinh ^ 



cosh f o^T^-^ 



sinh ( -2- 



M/ cosh 



sinh ( -2- 



(20) 



(21) 



(22) 



These are the solutions of the hnearized Eilenberger equa- 
tion in the ferromagnct. The corresponding expressions 
for are given by Eqs. (|Alip - (|A13p in the Appendix. 

Let us consider how these expressions reproduce the 
hmiting cases of strong and weak scattering, for which 
the solutions are well-knownii^. 



The boundary condition at x = read; 



,28 



Ao 



(25) 



Substituting the expression (|24p into the Usadel equation 
one can find the coefficient 



C. Dirty limit 

If the electron mean free path is the smallest charac- 
teristic length i.e. if <^ S,H-,Co,df, the frequent nonmag- 
netic scattering permits averaging over the trajectories^^. 
Then it is possible to write a closed system of equations 
for averaged functions. The linearization of the Usadel 
equations is allowed at the same conditions as for the 
Eilenberger equation. 

The linearized Usadel equation in the ferromagnct has 
the form 



{lu + iH) 



$+) = 



where ^? 



— £'//27rTc, and Df ~ ifVjjZ is the diffusion 
coefficient. The solution of Eq. (|23)) may be written as 



($+) = A cosh 



if 



iH 



■nTr 



(24) 



A = 



Ao 



(26) 



Let us consider the previously obtained averaged Eilen- 
berger function (jlSp . For small (.f the parameter q is 
large and 



1 + tf/i. + ltfKH 



q I TTm 

— arctan 



1 



TT^m^ + 







(23) It is seen that the sums in (ITil) and (|15p converge at 



2 (a; + iH) dj/df. At these values of : 



"^/32"^/ 



lijif < 1 



and 



VAT 



-^cosh 



iH 



7B sinh 



I 



-^(^coshf,:^ 



sinh ( — 



(27) 



where the suppression parameter describes the elec- 



tron transmission through the FS interface, see expres- 
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sions (|A14[) - (|A16p . After averaging over Of the expres- In the limit <C we have 
sion (g?!) coincides with ([Ml) with A given by It 
turns out that in the dirty hmit the first term of the ex- 
pression for the Eilenberger function plays the main 
role, as the second term reduces due to the spatial aver- 
aging (jisp . and the Eilenberger function = ($+). 



D. Clean limit 



— arctan ) w + ^^fUn 

Km \ Q 



For a larger electron mean free path £ / the sum in ()A4 
converges at 



jMf nm 1 

^ 7J ^ 1 — - — 
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and 




M/Cosh(^^^ MCOsh(2!Iifctf^i; 



sinh — 



sinh — 



sinh I — 



cosh 



(28) 



In this limit, the first term in the square brackets in 
Eq. p4p is small in comparison with the second one. 
Thus, one may conclude that the first term in Eq. (fT4)) 
plays the main role in the dirty limit, while the second 
term of Eq. (|14p describes mainly the clean limit when 
trajectories of motion are essential. 



The boundary conditions may essentially depend on the 
quality of the interface. In the case of a nonmagnetic 
specularly reflecting boundary between two metals these 
conditions read^ 

g"iRgf+g^_^)^Djglgl. (30) 



E. Josephson current 



9i = 92 = 5" 



(31) 



To calculate the Josephson current J of a SIFS ferro- 
magnetic tunnel junction we start from the expression^ 

J = 4eiV(0)w7rT^ ( {Imgn) cos9d{cos9) = 
= 4e7V(0)u7rT^ [ {lmg1^)iMln. (29) 



Here A'^(O) is the density of state at the Fermi surface, 
V is the Fermi velocity, gn is the matrix element of the 
antisymmetric part of the Eilenberger function that is 
define by the following expression: 



5" = - 
2 



9(0) -gH 
/+(^)-/+(- 



fio) - fi-o) 

~9{0)+9{-0) 



We apply the tunnel hamiltonian approach and use the 
boundary conditions on the dielectric interface. These 
boundary conditions are matching the quasiclassical elec- 
tron propagators g and / on both sides of the boundary. 



Here Di and i?/ = 1 — Dj are the interface transparency 
and reflectivity coefficients, the index 1(2) labels the 
functions on the right (left) side from the boundary plane, 
and the symmetric parts of the quasiclassical Eilenberger 
functions are defined by the equalities 

5^ = ^(5^-52^), 5^ = ^(5^ + 52^), (32) 



.91, 2(^') +5i,2(- 



-9i,2{S) - 9iA-9) 



Assuming that the insulating layer transparency is 
small, Dj ^ ~ 1, and taking into account that 

g1 ^ gf = 1, we can expand the boundary conditions 



(pO| in powers of Dj. 

For SIFS tunnel structures with s-wave pairing in the 
electrodes, in this limit it immediately follows that the 
functions g2{S) — gs{S) and f2{6) = fs{S) are indepen- 
dent on and coincide with the expressions for spatially 
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homogeneous superconducting electrodes 

ff = Ffie) exp {+iip/2} , f,=F, exp {~i^/2} , 

= , - gs = ; , . (33) 



^0 



A2 

^0 



Here Aq and ip are the absolute value and the phase 
difference of the order parameter in the electrodes. The 
boundary conditions pop reduce in this case to 



(34) 



^ £lf hf} - /2/1+ 2(.gi/2 - /i.g2) 
4 V2(/+g2-.gi/2+) ft h - fUi 

Then the expression for the Josephson current has 
the form 

'^'[/2+(^)/l(^) 

- /2(6')/i+(6')]cos6' d(cos6i). (35) 

Using the definition ([T|) and its symmetry properties ([TT]) 
and (jl2p , one can find the result for the Josephson current 
J of the tunnel junction J = Jq sunp, where 



J = -eN{0)v7TT ^ / 



(36) 

The Eilenbcrger function of the left electrode $+(/^s) is 
also defined by ([8]). 

The thin insulating layer is considered as a high poten- 
tial barrier for electrons. The transmission probability 
is inversely proportional to the exponent of a distance 
passed by an electron, i.e. ~ exp(— (i//cos0), where di is 
the dielectric thickness. Namely, it can be found from the 
well known expressions^ for the transmission coefficient 
of a square potential barrier of a thickness a, that 



Di{a) 



Akfkj + (fc2 + /c2)2 sinh^(fc2a) 



Akfki 



{kl + klY 



exp(— 2fc2a), 



(37) 



if k2a ^ 1. Here ki and ik2 are wave vectors of parti- 
cles outside and inside the barrier. The expressions for 
the transmission coefficient of the FS boundary ([T7 | . p^ 
were obtained in the same way. Taking into account that 
a = c?///i, one can write the dependence Di{^) in the 
form 



Dii^j.) = Do exp 



(38) 



Here a is a decaying coefficient, that depends on the 
thickness and material of the insulator. It is useful 
for a numerical calculation to redefine the value Dq = 
-Do exp (—a) , and finally 



Di{^) = Do exp -a 



1 



(39) 



III. DISCUSSION 



A. Main cases 



First, we start from the analysis of the Eilenberger 
function $+(a;,/Lt/), which describes the superconduct- 
ing properties of a ferromagnetic Josephson junction. In 
the simplest case ^+{x,fif) is an exponential function 
$+(x) oc exp(x/^), a; ^ 0, or a combination of exponen- 
tial functions, with a complex coherence lengtbii^ 



(40) 



Here ^1 describes the decay of the superconducting cor- 
relations at some distance from the FS boundary, while 
^2 defines the period of LOFF oscillation a^°'^^ . For the 
described FS structure, $+(a;, ///) is given by the expres- 
sions dHI) and ([20l)-(l22l)- Now we try to describe the 
main features of the proximity effect in an FS structure. 
If the superconducting layer is thick enough to be consid- 
ered as semi-infinite, and if in-plane nonuniformities con- 
nected with sample preparation and the domain structure 
of the ferromagnetic film can be neglected, the structure 
has one geometric parameter — the ferromagnet thick- 
ness df. The properties of the F material give three more 
characteristic lengths: the electron mean free path in 
the ferromagnet if, the nonmagnetic coherence length 
^0, and the characteristic magnetic length £^h- The main 
questions at this stage are: can the Eilenberger function 
be approximated by an exponential function, and what is 
the relation between Ci, ^2 and the characteristic lengths 
mentioned above. 

The Fermi velocity for usual ferromagnetic metals^i is 
Vf ^ 3 ■ lO^cm/s, and if even the temperature T ^ Tc, 
^0 ~ 100 nm, i.e. ^0 is much larger than other parameters 
of the problem in usual cases. The value of ^0 is diffi- 
cult to decrease. The gradual increase of ^/ leads to the 
following four cases: 

1. if <^ (,H,df — the well investigated dirty limit. 
This condition allows averaging of the Eilenberger func- 
tions over trajectories and using the Usadel equations. 

is given by the expression (p4)) . and for a strong ex- 
change field {H > fcsTe), 6=6 = y/2if^H/3. 

2. if ~ ^H,df — the intermediate case. The Us- 
adel equations cannot be used, and one needs to solve 
the Eilenberger equations taking into account the $+(6') 
dependence. Up to now the linearized equations were 
solved only for an SFS junction^i^, and the obtained ex- 
pressions are so complicated that they are very difficult 
to analyze. This is the most interesting case for our anal- 
ysis, which allows to find out when this case reduces to 
the dirty limit and under which conditions the Usadel 
equations are applicable. 

3. < df,£,o — the clean limit. Here, $+(a;,0) 
is defined mainly by the second term of Eq. ([T4| , and^i^^ 

4. if df — the ballistic regime cannot be considered 
within the framework of our approach, excluding the case 
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^1 <S^ df when multiple reflections from the interfaces 
can also be neglected. The ballistic SFS junction was 
considered earlier— i^^'^^''^^" —, and the dependence Jc{df) 
cannot be presented as an exponent in the general caso^^. 



berger function 1), and in further calculations we 

will also use the form D{^f) = £'2/^/ for the transmission 
coefScient of the FS interface, as it is most applicable to 
materials used in experiments. 



B. Analysis of cases 1-3 

In all plots presented below we use the same normaliza- 
tion for all lengths: x, df, If, ^0, £,h, £.2- All of them 
are normalized to some unit length S, which, in fact, can 
be chosen arbitrarily, e.g. S = 1 nm. Note that only the 
ratios between the lengths given above are important, 
i.e., if one divides the above list of lengths by any arbi- 
trary constant S all the results remain unchanged. For 
example, in all plots we use ^0 =81. 

a. Influence of different types of FS boundary trans- 
parency. Figure [T] shows the spatial distribution of the 
functions calculated from Q for = 7, £h = 3, 

df = 10 and Jc{df) calculated from ([55)1 for = 
7, = 3, and the insulating layer decay parameter 
a = 5. It is clearly seen that the angular dependence of 
SF interface transparency influences both, the amplitude 
and period of oscillations of the presented curves. For 
the particular choice of parameters this influence can be 
easily explained. 

For the case of a transparency which is independent 
on fj,f, the Eilenberger functions are initiated from the 
SF interface in all directions, thus leading to the largest 
amplitude of oscillations. Simultaneously, the 

contributions to ($+(a;)) coming from rapidly decaying 
Eilenberger functions with a large argument 9f result in 
the appearance of the smallest period of oscilla- 
tions. The stronger the angular dependence of the SF in- 
terface transparency coefficient, the smaller are the con- 
tributions to average values from the rapidly decaying 
Eilenberger functions, and as the result the larger is the 
period and the smaller is the amplitude of ($+(.t)) os- 
cillations. It is necessary also to mention that ($+(.t)) 
decays more rapidly with |x| than 1). This differ- 

ence is the smaller the smaller is // . 

From (j36|) it follows that in general the contribution 
to the junction critical current comes not only from 
$+(a;, 1), but also from <I>+(x,/i/) located in a narrow 
domain of /i/ nearby /i/ = 1. Due to that, the Jc{df) 
curve appears to be sensitive against the angular depen- 
dence of the SF interface transparency. We see that the 
stronger this dependence the smaller is the amplitude of 
oscillations and the larger are the distances of the to tt 
transition points in Jc{df) from the SF interface, while 
the period of Jc{df) oscillations is practically insensi- 
tive against the form of D{fif) and coincides with that 
of 3>+(a;,l). If the ferromagnetic layer is thick enough, 
the sum over uj in ([55]) converges rapidly, and the main 
contribution is given by the first term of the sum, which 
is determined by the real part of the Eilenberger function 
1). For these reasons we will focus below on the 
examination of the properties of the real part of the Eilcn- 
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FIG. 1: (Color online) (a) Averaged Eilenberger function 
{^+{x)) and (b) Jc{df)RM product for different types of 
the FS interface (|16P - (|19p . The F-layer parameters are If = 
7, = 3 (a,b), df = 10 (a), the insulating layer decay pa- 
rameter is a = 5 (b); T ~ 0.5Tc. 

b. Spatial profile of the Eilenberger function. The 
direct comparison (see Fig. [5]) of the results of numerical 
calculations of the real part of the Eilenberger function 
^+{x, 1), making use of Eq. (|2ip . with the curves which 
follow from the simple analytical expression 



(x,l) 



cosh(d//C) 



cosh ■ 



■df 



(41) 



confirms that it is possible to approximate $+(a;,l) by 
this simple formula. It satisfies the boundary condition 
on the dielectric interface ^ and contains the factor 
$+(0,1), which has been calculated numerically from 

The quality of the approximation (|4ip was checked for 
different values of parameters for all dependencies inves- 
tigated further. It was found to be satisfactory in the 
whole accessible range of parameters df > if 01 £1. This 
approximation yields the values of ^, i.e. ^1 and ^2- 

c. Dependence of £ on the junction parameters. The 
attempt to find ^ in the area of intermediate scattering 
was made earlier—. It was supposed in Ref. [s^ that 
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FIG. 2: (Color online) The exact ^+{x,l) (solid lines) and 
approximated ^+{x, 1)°pp (dashed lines) Eilenberger function 
(a) for different values of corresponding to a local maxi- 
mum, minimum, and smooth region of the dependence £,i{^h)', 
see Fig[3l^b), at £/ = 7, d/ = 10; (b) for different thicknesses 
df at £f = 7, £,H = 3.8. T f» 0.5Tc. 
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may be found solution of the equation 

= arctan(i£//A:/^), (42) 

which determines the poles of the sum in Egs. p^ . ([15]) 
for ^+{x) and ($+(x)), respectively. Unfortunately, the 
inverse tangent is a multivalent function. In the general 
case this fact prevents to represent ^+{x) and ($+(0;)) as 
only the sum of residues of summable functions in (|14p , 
(jlSp . An exception occurs in the range of parameters 
for which the ferromagnet is close to the dirty limit. In 
this case the sum in (|14p and (|15p converts faster than 
the multivalent nature of the inverse tangent becomes es- 
sential. These difficulties had been first pointed out in 
Ref. [131, where it was also demonstrated that the solu- 
tion of the Eilenberger equations used in Ref. [s^ does 
not transfer correctly to the solution of the Eilenberger 
equations in the clean limit. 

To avoid these mathematical difficulties we have de- 
cided to find out the dependence £,(£,H,if) from the direct 
fitting of the real part of ^+{x) calculated numerically 
from ^4]) by approximating formula (|4T|). The result of 
this procedure is demonstrated in Fig. [31 It shows the 
dependencies of ^1 and ^2 on the magnetic length at 
different £f. The solutions from Ref. [s^ are also pre- 



FIG. 3: (Color online) Decay length ^1 and oscillation param- 
eter ^2 vs magnetic length for different values of and 
df, that are denoted in the plots. The corresponding values 
of ^ from Ref. [3^ are also shown and marked by *. 

sented in Fig. [3] for comparison. They coincide with our 
values wheir ^ if i-e. in the dirty limit. We have 
taken two cases for the numerical investigation: if = 1, 
that is close to the dirty limit, and if = 7, that is close 
to the clean case. 

It is seen, that ^1 oscillates with the exchange energy 
H (so as £,H ~ an interesting result, that was not 

noted earlier. These oscillations appear in rather clean 
samples if > df/n and for < df/n. They origi- 
nate mainly from the second term of the expression (jl4p . 
which plays the main role in the clean case. The func- 
tion £,{£,h) changes with increasing df (cf. Fig. [3^b-d)). 
At large enough df, electrons have a time to scatter in 
the ferromagnet, their trajectories shuffie, and the dom- 
inant role of the "clean" term vanishes, the oscillations 
disappear, see Fig. ^d). We note that the smaller is 
the probability of dcphasing of the Eilenberger functions 
due to electron scattering in the F layer, the stronger are 
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the interference effects between the IF and FS interfaces. 
This interference demands, for given values of and 
the fulfillment of the boundary conditions at these in- 
terfaces, which fix the value of the phase derivatives at 
X = and x = df. As it follows from Eq. ([20|)-(l22]), these 
conditions cannot be satisfied on the class of monotonic 
^1 and ^2 functions. Therefore the transition from the 
dirty to the clean limit should be accompanied by a non- 
monotonic behavior of and ^2 as a function of for 
fixed £, see Fig. [31 or as a function of £ for fixed ^h, see 
Fig.H 

In the dirty limit the spatial oscillation length ^2 ~ 
^y2£f£Jf/3, while in the clean limit ^2 ~ Near the 
point marked by X in Fig. [3] there is a crossover between 
the dependencies corresponding to dirty and clean lim- 
its, respectively, i.e. this region can be considered as a 
boundary between the clean and the dirty cases. 
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FIG. 4: (Color online) Decay length ^1 (solid lines) and oscil- 
lation parameter ^2 (dashed lines) vs electron mean free path 
if in the ferromagnet for different values of ^h; here df = 10, 
T«0.5Tc. 

The dependence of ^ on the mean free path £ for differ- 
ent £,H is presented in Fig. Usually the decay length 
^1 as well as the oscillation period 27r^2 increase with 
the mean free path, which corresponds to the results in 
Ref. [H. However, the dependence £,{£f) may be non- 
monotonic for some values of the ferromagnet thickness. 

In the presented figures one can see, that ^1 may be 
larger or smaller than ^2 and may behave nonmonotoni- 
cally. This depends not only on the material constants of 
ferromagnetic material, but also on the thickness of the 
F layer. 

The found values ^1 and ^2 together with the simple 
dependence (|4ip could be widely used for various estima- 
tions, approximate calculations, and fitting experimental 
data, for example for a measurement of the density of 
states (DOS) in an FS bilayer or calculation of thick- 
ness dependence of the critical current of SIFS Josephson 
junctions. 

d. Jc{df) in the limit of large df. As an example, 
we consider the critical current of SIFS Josephson junc- 
tions in the limit of large F layer thickness df ^ In 



this limit the main contribution to the critical current in 
is given by the term at w = ttT. Suppose further that 
only electrons, that are incident in the direction perpen- 
dicular to the FI interface, provide the current across it, 
from dSe]), dHI) we have 



Jc = — ^ — Re 

eRN 



Ao$+(0,1) 



v/(7rT)2 + A-^cosh(d//C) 



(43) 



The comparison with the exact numerical calculation 
has shown that this expression ([43]) well approximates 
Jc{df) starting from df > ^1. The thickness depen- 
dence of Jc calculated from at the value ^ taken at 
df = 50, see Figl3];d), for If = 7, = 2.1, T « O.bTc 
is shown in Fig. [5] by the solid line. The dashed and 
dashed dotted curves in Fig. Ogive the results obtained 
by numerical calculation, which have been done for the 
same parameters with the use of the exact expression 
for $+(— d/, /i/) in ((36|) and at different insulating layer 
thicknesses described by the parameter a. 

It is clearly seen that the larger is a the closer is the 
the result to the approximation formula (|^ . We may 
also conclude that the difference between the asymptotic 
solid curve and the dashed curves in Fig. [5] calculated for 
finite values of a occurs only in amplitude and positions 
of the to TT transition points, while the decay length and 
period of Jc{df) oscillations are nearly the same. This 
means, that if in experiment we are mainly interested in 
the estimation of the ferromagnet material constant, ^, 
we may really use for the data interpretation the simple 
expression and consider the coefficient $+(0,1) in 
it as a phenomenological parameter, which depends not 
only on the properties of ferromagnetic material, but also 
on in a rather complicated way, see Figl^Jb), as well as 
on the form of the transparency coefficients at the FI and 
FS interfaces. Although our results for SIFS junctions 
give the same ^1 and ^2 as for SFS junctions (at least 
for a thick F layer), different boundary conditions result 
in different order parameter amplitude and phase. This 
results in a different Jc{df) dependence, similar to the 
results obtained earlier in dirty limit^. 

It is interesting to note that in the case of a rather 
thin insulating barrier and rather clean ferromagnet, 
^1 and ^2 as measured from the critical current of SIFS 
junctions may differ from ^ as measured from the DOS 
on the free F surface of an analogous FS bilayer. This is 
because the DOS measured by low-temperature scanning 
tunneling spectroscopy (which may only yield contribu- 
tions from electrons with normal incidence) is directly 
given by $+(— d/, = 1), while Jc includes contribu- 
tions from other angles. 

e. Applicability of the Usadel equation. The condi- 
tion that allows using the Usadel equation is well known 
— this is strong nonmagnetic scatteringi^, namely, ^/ <C 
£.H,df,^n- What does the symbol mean exactly? 

In order to illustrate this, Fig. [6] shows the dependence 
Jci^f) given by the expression (pij) . Two parts of (PT|) . 
called for convenience " dirty" and " clean" , are given by 
the first and the second terms of the expression (PT|) . The 
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FIG. 5: (Color online) JcRn product as a function of the 
ferromagnet thickness df at different insulating layer thick- 
nesses described by the parameter a, and the approximation 
by expressions (|36|l . (|41| l at the value ^ taken at d/ = 50, see 
FigHd), here If = 7,^h = 2.1, T 0.5T,,. 
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FIG. 6: (Color online) JcRn product vs ferromagnet mean 
free path If {£_h = 10, T O.STc, a = 5) for df = 10 (a) 
and df = 30 (b). Each plot shows JcRn obtained from the 
Usadel function pSI) . from the Eilenberger function (|21|) . 
as well as from only the first "dirty" term or only the second 
"clean" term of (|2T1). 



critical current density as given by the Usadel function 
([27|) . describing the dirty limit, is also shown. Jc(^/) 
coincides with the Usadel solution at ^/ « 0.1£^h, then 
the first term in square brackets in Eq. (|2ip dominates 
and the /ij-dependence becomes nonsignificant due to 
the spatial averaging as a result of multiple scattering. 
This means, that the Usadel equations become appropri- 
ate to use if the parameter Hrf ^0.1. This result was 



intuitively clear, but demanded a proof due to plenty of 
investigations, where the Usadel equations was used at 
HTf < 1. 



IV. COMPARISON WITH EXPERIMENT 

There are only few experiments on SIFS tunnel junc- 
tions with moderate scattering in the F-layer-i2i2S. In 
Ref. [lo| experimental data points for Jc{df) are rather 
sparse, while in Ref. [39j the density of data points per 
period is much larger. We thus decided to compare our 
theory to the data of Ref. [s^ . 

This article^ also contains attempts to fit the ex- 
perimental data by theoretical curves, however, without 
much success. Some of these attempts, see Fig. 5(a) in 
Ref. [1^, use a theory developed for a very clean (bal- 
listic) SFS junction. The predicted oscillation period of 
Jc{df) was significantly smaller than the one in the ex- 
periment. Figure 5(b) in Ref. [s^ shows the same ex- 
perimental data together with fits using dirty limit the- 
ory. The best fit was achieved with = 0.66 nm and 
^2 = 0.53 nm. However, the theory for a dirty ferromag- 
netic junction cannot explain > ^2, and all of these 
theories, both clean and dirty, do not take into account 
the insulating layer. 

In Ref. [33 and Ref. [i^ it was found, that the mag- 
netic anisotropy of the F layer changes from perpendicu- 
lar to in-plane with an increase of the F-layer thickness. 
Perpendicular magnetic anisotropy in a polycrystalline 
thin film may occur due to several reasons: the mecha- 
nism described by Neel^ related to the absence of near- 
est neighbor F-layer atoms near interfaces, the magneto- 
striction mechanism, the surface roughness, and the one 
associated with microscopic-shape anisotropj*^. The lat- 
ter may also vary for samples sputtered under different 
angles^. Perpendicular magnetic anisotropy can change 
to an in-plane anisotropy due to a competition with the 
shape anisotropy of a film. However, none of the mech- 
anisms listed above can be expected to give a significant 
jump of the exchange magnetic energy H at the transi- 
tion between the two types of anisotropy. There is some 
shift of the experimental curve, see Fig. 2(a) in Ref. [soj . 
at df = 3.8 nm, where the anisotropy changes. But there 
are no reasons to assume, that this change of anisotropy 
is associated with a significant change of H . Therefore, 
we try to describe the whole experimental curve Jc{df) 
using the same value of H . 

The common problem for all explanations of this ex- 
periment is the location of the first minimum of Jc{df) 
at a rather large df ~ 3 nm in comparison with the oscil- 
lation length ^2 ~ 0.5 nm. The authors^ attribute this 
to the existence of a nonmagnetic dead layer with rather 
large thickness ddead = 2.26 nm in the Ni layer. 

We note, that ddead, in general, depends on df. We 
assume that ddead = df up to some value dmax and 
then either ddead stays constant at a saturation value 
^dead ~ d^nax with further increase of df, or it is reduced 
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to a constant value rf^"*^ < dmax- For the latter case 
there are the following reasons. The dead layer may con- 
sist of a nonmagnetic Ni alloy with magnetic clusters or 
islands, which are not connected by the exchange inter- 
action. A very thin film of a ferromagnetic material may 
be paramagnetic. When one covers it by extra magnetic 
layers, some part of its thickness may magnetize. First, 
magnetic clusters uncoupled by the exchange interaction 
can couple if they are covered by extra magnetic layers. 
Second, a thin homogeneous film of a magnetic material 
of the thickness of a few monolayers exhibits nonmagnetic 
properties alone, but becomes magnetic entirely when its 
thickness increases, due to the transition from the two- 
dimensional to the three-dimensional case. 

The dead layer is described as a normal layer, that 
yields Jc{df) ^ exp( — |^). The best fitting presented 
in Fig. [71[a),(b) by straight short-dashed lines, gives the 
valued = 0.68 nm. The normal metal coherence 
length S.N depends both on the nonmagnetic and para- 
magnetic scattering. The latter may be significant in 
the dead layer in the presence of many magnetic inho- 
mogcneities (clusters and impurities). Following a the- 
ory developed for SFS junctions^ taking into account 
magnetic scattering, one can find the expression for the 
coherence length of a normal metal (at H = 0) with 

magnetic impurities ~ \J 2(TrT+i/T } ■ '^here r,„ is the 
magnetic scattering time. We assume that the mean free 
path has the same value as in the rest of the magnetic 
Ni layer, £f = 0.6 nm, that is given by the fit. This 
yields t„i = 1.68 • lO^^^'s. For comparison, the nonmag- 
netic scattering time is found to be almost 8 times less. 
Fig. [7Ka),(b) shows that the obtained value of is very 
close to the decay length of the magnetic SIFS junction 
at df > dmax- In the absence of the exchange field its 
pair-breaking role is played by the magnetic scattering in 
this structure. 

We consider two possibilities to explain the experimen- 
tal datai^. First, we assume that the detected minimum 
of Jc{df) is the first one, see Fig. [7|a). Treating the 
dead layer as a normal layer (H = 0), our best fit yields 
d'fj^^ = dmax = 2.37 nm. The fitting yields the fol- 
lowing values of the F-layer parameters: £f = 0.6 nm 
and = 0.6 nm, that corresponds to H = 1782 K or 
H = 154 meV. 

As a second scenario we assume that the detected min- 
imum of Jc{df) is the second one, see Fig. [TJb). The 
first Jc{df) minimum may not show up if the effective 
thickness of the dead layer decreases as soon as df is 
above the threshold value dmax- Then, at the value oi df 
where Jc{df) is supposed to have the first minimum, the 
F layer may still be nonmagnetic, while at large values 
oi df it has some small value d^e^d- For that scenario 
'^Tead = 0.2 nm provides the best fit. The mean free path 
was found as £f = 0.6 nm, = 0.956 nm, that corre- 
sponds to H = 1119 K or H = 96.4 mcV. 

In both fits we used the following parameters. The 
temperature T = 4.2 K and Tc = 8K were taken like 



in the experiment^, that yields = 81 nm. Fermi 
velocities^! Vf = 2.8 • W cm/s and^ Vs = 6.2 ■ lO'^cm/s 
were taken. Since Vs > Vf, we use the FS boundary con- 
dition in the form ([T7)) . The thickness of the insulating 
AI2O3 barrier was not measured exactly; still the value 
a = 5 seems to be realistic, see Eg. ([5^1). The values of 
the exchange field for both cases are within the range 
obtained for Ni in SFS experiments^ '^^i^^-^^ . The rela- 
tion £f ^ does not allow using the Usadel equations, 
however, this is not far from the dirty limit. 

Comparing both fits shown in Fig. [7Ka),(b) we cannot 
give a definite answer which minimum (the first or the 
second) was observed in the experiment^, but we are 
able to explain the fact that ^1 > ^2 in any case. Such 
relation between ^1 and ^2 cannot be obtained from a 
dirty limit theory. 

The Josephson junctions used in the experiment^ con- 
tain an extra normal layer of Cu between the I and F 
layer. A normal layer cannot change ^, which mainly 
depends on properties of the magnetic layer. However, 
it may change the boundary conditions. Therefore, our 
theory cannot be directly applied for the determination 
of the exact value of the dead layer, but may explain the 
obtained value of ^. 

For the analysis of the first minimum position and the 
dead layer estimation we may also use the theory devel- 
oped earlier for the dirty limit, since the estimated value 
of ^2 • From the analysis of the uniform part of the 
solutions obtained in Ref. gl] for dirty SINFS and SIFNS 
junctions, it follows that the first minimum of the Jc{df) 
dependence must be at [df — ddead)/S,2 ~ 37r/4. If we as- 
sume that the detected minimum of Jc{df) is the first 
one, then, the estimation yields = dmax = 1.8 nm. 

This value is rather large (it is larger than the oscilla- 
tions halfperiod), similar to Ref. [33! ■ 1^ assume that 
the detected minimum of Jc{df) is the second one, then 
the first one must be located at df ~ 1.4 nm and the 
dead layer thickness d^"*^ = 0.2 nm. There is a small 
shift of experimental data for Jc{df) at the thickness 
d/ « 1.4 nm, but a real O-tt transition was not detected, 
possibly, due to the fact that d^g*^; < dmax- Generally 
speaking, if we assume that ddcad{df) jumps from dmax 
to djcad c'max then ouc should see a jump on Ic{df) 
dependence. However in our case it not observed, see 
FigEJa) and (b). We suppose that this is related to the 
fact that « ^1. 

Figure [TJc) contains also calculated curves, plotted for 
the same model parameters as in Fig. [TJ^b) but for a 
cleaner ferromagnet. Usually the purpose of the exper- 
imental investigations of SIFS structures is the design 
of a TT- Josephson junction with a large JcRn- Rn is 
determined mainly by the insulating barrier. Our calcu- 
lations show that by increasing I f one can increase Jc by 
1-2 orders of magnitude, see FiglTl^c). It would also be 
reasonable to delete extra normal layers to achieve the 
TT-phase at a smaller df, and consequently, to have larger 
Jc; see also Refs. fsslli^ . 
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FIG. 7: (Color online) (a) and (b) show a comparison of 
Jc{df) from the presented theory with experimenli^. Dashed 
lines are fits to the experimental data in the range d/ < 
2.37 nm, where the F layer can be interpreted as nonmagnetic, 
(c) Dependencies Jc(d/) at lower scattering. Parameters are 
like in (b). The regime for which multiple reflections cannot 
be neglected are marked by dotted lines. 



ness df, and can be nonmonotonic as a function of or 

The approximation of Jc{df) by an exponential fmic- 
tion or its combinations has some restrictions. The ap- 
phcabihty of the Usadel equation has been established 
for £f < O.l^H- 

The developed approach has been used to fit exper- 
imental data^ providing a satisfactory fitting of the 
Jc{df) dependence. It allows to explain the values 
Ci > ^2 and to give some practical recommendations on 
how to increase the JcRn product. 
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Appendix A 



Substituting the expression (|13p into the Eilenberger 
equation ([5]), multiplying the obtained equations by 
cos [Trkix + df)/df] , and integrating them over x one can 
find the relation between coefficients Qm and B (|10p . 
namely 



{Qr, 



B 



df M 



sinh 



Cf 



kfM 



(Al) 



where 



M = 7rW^ + l, /.^ 



q mn 
— arctan 



\ M I -nm q 
Averaging both sides of ()Aip over the angle we get 



V. CONCLUSION 

The SIFS ferromagnetic tunnel Josephson junction has 
been investigated within the framework of the quasiclas- 
sical Eilenberger equations, that allow a description of 
both, the clean and the dirty limits, as well as an ar- 
bitrary scattering. The Eilenberger function $_|_(a;,/^/) 
may be approximated by the simple formula (j4ip within 
the entire range of considered parameters. The decay 
length ^1 and the oscillation period 27r^2 depends not 
only on the mean free path and the exchange energy 
H in the ferromagnet, but also on the ferromagnet thick- 



(0 )-( B^iziTsinh^ 



1 

kfM 



1 

kfM 



(A2) 



A substitution of (jA2[) into (jAl[) finally gives the relation 
between the coefficients in expression ([Tc 



Qr. 



B 



df M smn 



M 



[kf-{^)] ■ 



(A3) 



Expressions (jA3p and ([T^ permit to rewrite the solution 
of the Eilenberger equations in the closed form 



13 



E 



■ COS 



'!Tm{x + df) 



m— — cjo 



MiH) [kf-iii)] 
and for an isotropic component ($+) of this solution to get 

, (-!)'" M 
\ -"^P-j 



^(M)T7rx-sinh^ 



B(fif) cosh ( q 



■Km{x + c?/) 



(A4) 



(A5) 



The form of the solution of the Eilenberger equations in the S layer essentially depends on transport properties of 
the FS interface. 

If the transparency D{fif) is small, then in the first approximation on D{^f), we can neglect the suppression 
of superconductivity in the S region and consider the order parameter and the Eilenberger functions as constants 
independent on space coordinates, which are equal to their bulk values, thus leading to the boundary condition for 
determination of integration constants -B(/x/) in the form of Eq. ((TU]). From the boundary condition pU)) we get 



An 



1 



(A6) 



and obtain the analytical solution in the form 
analysis than solutions previously used in Refs 

For the described three forms ([T5)) - p^ of D{iif) dependencies we get 



VA^+w2sinh(^) ' 
that, from our point of view, is more convenient for the further 
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1 arctan 
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(A8) 

(A9) 
(AlO) 



Substituting the averages (|A7[) - (|A9[) into the expression we get the Eilenberger function for different FS interface 
transparencies The functions (pOjl -ip ^ averaged over the angle have the following forms. 

First, for D{^if) = Di 



2q [kf 



f - 1) ^ ^ kf-r^ arctan ^ 

1—1-' 



COS ■ 



df 



(All) 



Second, for D{iif) = /iZ?2 

Third, ior D{nf)^ fi^Ds 
(*+) 
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Our calculation yields the following value of the sup- £f {^D{^)) ^ /i^J df /2'kTc. It differs from the one ob- 
pression parameter for the Usadel equation: 7^ = 
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tained earlj«2^ by a factor 2. We suppose, that this is 
a consequence of an approximation we have made in 
boundary conditions (fTO|) at the FS interface. We have 
neglected all spatial variations in the S part as well as the 
back influence of the ferromagnet on the superconductor, 
and this factor 2 is the price for this approximation. If 
wc use the boundary condition (1101) in its full foim^, we 
get for jB- 

7b = ?^^(m^(/^))-'. 



27rTc 

Our three models of the FS interface yield 

2 ^/ /,A-i_4 if 
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if^y'^T, — -j==, (A14) 



2itTc ^-^V 27rT, 
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3 f^^'^')"=2^f^' 
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(M^)-^ = |^f=. (A16) 
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